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Abstract. We propose and study discontinuous Galerkin methods for strongly 

degenerate convection-diffusion equations perturbed by a fractional diffusion 

(_P ' (Levy) operator. We prove various stability estimates along with convergence 

Ch ' results toward properly defined (entropy) solutions of linear and nonlinear 

pH , equations. Finally, the qualitative behavior of solutions of such equations are 

j^ ■ illustrated through numerical experiments. 
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1. Introduction 

We consider degenerate convection-diffusion equations perturbed by a fractional 
diffusion (Levy) operator; more precisely, problems of the form 



< 

^' : \ut + f{u), = {a{u)u,),+bC[u] (x,t)eQT = Kx (0,T), 

C^ [ \u{x,0) = uo{x) X e R, 

where /, a : M — ?► R (a > and bounded) are Lipschitz continuous functions, b > 
is a constant, and £ is a nonlocal operator whose singular integral representation 
p^ ; reads (cf. EZIIE]) 

t^ ; C[u{x,t)]^cx — . '^ ^ ' ' dz, Xe (0,1) andcA > 0. 

"^ ' For sake of simplicity, we assume /(O) = 0. The initial datum uq : K — >■ R is 

chosen in different spaces (cf. Theorems 14.21 14.41 and I5.8p depending on whether 

|;0 I the equations are linear or nonlinear. 

The operator C is known as the fractional Laplacian (a nonlocal generalization 
of the Laplace operator) and can also be defined in terms of its Fourier transform: 



(1.2) c[u{;tm)^~m'uiu). 

As pointed out in [2j[T2j[27], u{-,t) has to be rather smooth with suitable growth 
$H ' at infinity for the quantity C[u] to be pointwise well defined. However, smooth 

solutions of (II. ip do not exist in general (shocks may develop), and weak entropy 
solutions have to be considered, cf. Definition 15.11 and Lemma [A. II below. 

Nonlocal equations like (jl.ip appear in different areas of research. For instance, 
in mathematical finance, option pricing models based on jump processes (cf. [8]) 
give rise to linear partial differential equations with nonlocal terms. Nonlinear 
equations appear in dislocation dynamics, hydrodynamics and molecular biology 
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[13]; applications to semiconductors devices and explosives can also be found p5] . 
For more information about the possible applications of such equations we refer the 
reader to the detailed discussions in [2], |3], and [TT| . 

Equation (jl.ip consists of three different terms: nonlinear convection f[u)x, 
nonlinear diffusion {a{u)ux)x, and fractional diffusion C[u\. It is expected that the 
effect of a diffusion operator is that solutions become smoother than the prescribed 
initial data. In our case, however, a can be strongly degenerate (i.e., vanish on 
intervals of positive length), and hence solutions can exhibit shocks. We refer to 
[14l[T3] for the case when 5 = 0, and to 0[5j for the case when A € (0, 1) and a = 0. 
The issue at stake here is that the fractional diffusion operator may not be strong 
enough to prevent solutions of (jl.ll) from developing discontinuities. However, and 
as expected, in the linear case (/(m) = cu, a{u) = au with c G M, a > 0), some 
regularity can be proved (cf. Lemma l4.ip . 

An ample literature is available on numerical methods for computing entropy 
solutions of degenerate convection-diffusion equations, cf. [3 [T31 HH HSl HZl [13 [2S1 
[231 [10]. To the best of our knowledge, there are no works on nonlocal versions 
of these equations. However, for the special case of fractional conservation laws 
(a = 0) there are a few recent works [101 [TTJ [S] . Dedner and Rohde [TU] introduced 
a general class of difference methods for equations appearing in radiative hydro- 
dynamics. Droniou [11] devised a classs difference method for (II. ip (a = 0) and 
proved convergence. Cifani et al. [F applied the discontinuous Galerkin method to 
p.ip (a = 0) and proved error estimates. Finally, let us mention that the discontin- 
uous Galerkin method has also been used to numerically solve nonlinear convection 
problems appended with possibly nonlocal dissipative terms in [^ [^^ . 

The discontinuous Galerkin (DG hereafter) method is a well established method 
for approximating solutions of convection [6 and convection-diffusion equations 
[7[ [5^ . To obtain a DG approximation of a nonlinear equation, one has to pass 
to the weak formulation, do integration by parts, and replace the nonlinearities 
with suitable numerical fluxes (fluxes which enforce numerical stability and conver- 
gence). Available DG methods for convection-diffusion equations are the local DG 
(LDG hereafter) 7 and the direct DG (DDG hereafter) 20 . In the LDG method, 
the convection-diffusion equation is rewritten as a first order system and then ap- 
proximated by the DG method for conservation laws. In the DDG method, the 
DG method is applied directly to the convection-diffusion equation after a suitable 
numerical flux has been derived for the diffusion term. 

This paper is a continuation of our previous work on DG methods for fractional 
conservation laws 5 . We devise and study DDG and LDG approximations of (|l.ip . 
we prove that both approximations are L^-stable and, whenever linear equations are 
considered, high-order accurate. In the nonlinear case, we work with an entropy for- 
mulation for (jl.ip which generalizes the one in [301 [13]; and we show that the DDG 
method converges toward an entropy solution when piecewise constant elements are 
used. To do so, we extend the results in [T3] to our nonlocal setting. Finally, we 
present numerical experiments shedding some light on the qualitative behavior of 
solutions of fractional, strongly degenerate convection-diffusion equations. 



2. A SEMI-DISCRETE METHOD 

Let us choose a spatial grid Xi = zAx (Ax > 0, i G Z), and label Ii = {xi^Xi^i). 
We denote by P*'{Ii) the space of all polynomials of degree at most k with support 
on /i, and let 



V'^^iv: v\i^ e P'^ih), % e Z}. 
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Let US introduce the Legendre polynomials {ipo,i, 'Pi.t, ■ • • j Vk,i}, where (pj^i G P^{Ii). 

Each function in P^(Ii) can be written as a linear combination of these polynomials. 

We recall the following well known properties of the Legendre polynomials: for 

all i € Z, 

r ( ^^ for p = q 

j^ fp^^fq,^dx = | 29+1 ^^^^^^.^^ , fpAx^+i) = 1 and ipp4x+) = (-l)P, 

where (p{xf) = lim^^^± ip{s). 

The following fractional Sobolev space is also needed in what follows (see, e.g., 
[1] or [H Section 6]): 

II l|2 _ II ||2 I I |2 

ll"ll_H"^/2(R) — II"IIl2(R) + P|^A/2(R), 

with semi-norm IuI^a/2/r-i = Im Im \z-x\^+>' dzdx. Finally, let us introduce the 
operators 



[p{xi)] ^p{xj)-p{x.^ ), p{xi) = -{p{xj) +p{x^ )). 

From now on we split our exposition into two parts, one dedicated to the DDG 
method and another one dedicated to the LDG method. 

2.1. DDG method. Let us multiply (|l.ip by an arbitrary v E P'^{Ii), integrate 
over li, and use integration by parts, to arrive at 

(2.1) '^- ^^' 

+ / a{u)uj:Vx-h(ui+i,Ux^^+i)v^^-^+h{ui,u^^i)vf = b j C[u]v, 

where f{ui) = f{u{x{)), h{u,Ux) = a{u)ux and {ui,Ux^i) = {u(x{),Ux{xi)). Let us 
introduce the Lipschitz continuous E-flux (a consistent and monotone flux), 

(2.2) f{u,)^f{u{xr),u{xt)). 

Note that since / is consistent (/(u, u) — f{u)) and monotone (increasing w.r.t. its 
first variable and decreasing w.r.t its second variable), 

(2.3) / [!{x)-f{u-,ut)\dx>Q. 

J u. 

Following Jue and Liu ^D], let us also introduce the flux 

h{ui) = h{u{x-), . . . , d^u{xl),u{xj), . . .,d^u{x+)) 

r./ XI Lfc/2J 

= /3o^^^+A(^.)x+ E PmAx'"'~'[dl"'A{u,)], 
Ax ^-^ 

m—l 

where A{u) = J a and the weights {/3o, . . . , P[k/2\ } fulfill the following admissibility 
condition: there exist 76 (0,1) and a > such that 

(2.4) ^/,(„,)[^,]>«^I^K]-7^ /'a(«)(z.,)2. 

Note that the numerical flux h is an approximation of A{ui)x = a(u{xi))ux{xi) 
involving the average A{ui)x and the jumps of even order derivatives of A{ui) up 
to TO = k/2. For example, if fc = and (3q — 1, then 

M...)-^[^(^.)]- ^^"^"-'\~/^"^"-^\ 
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and this function satisfies condition (|2.4I) . In this case (fc = 0), 



A{ut)x = a{u{x,))d^u{xi) = -ia{u{x:l))dxu{xf) + a{u{x. ))dxu{x^ )j = 0. 
When k > 2, some extra differentiabihty on a is required. For example, with k — 2, 

Lfc/2J 

J2 PmAx^"'-'[dl'''Aiu,)] = (3,Ax[dlA{u,)] = l3iAx[a'{u,){d,u,f + aMd^u,]. 

m—l 

We see that the flux h is locally Lipschitz if a is sufficently regular, and that h{0) = 
for all k. Let us rewrite (|2.ip as 



utv - I fiu)vx + fiui+i)v-^^ - f{u{)vf 

+ / a{u)uxVx - h{ui+i)v^j^^ + h{ui)vl =b I £[u]v, 



and use the initial condition 

(2.6) / u{xtQ)v{x) dx — / UQ{x)v{x)dx. 

The DDG method consists of finding functions u : Qt — > K, u{-,t) e V'' , and 

fc 

(2.7) u{x, t)=Y,Yl UpAt)VpA^), 

which satisfy ^I^-^X^ for all v £ P''{1,), i £ Z. 



2.2. LDG method. Let us write a(u)Ux — y^a{u)g(u)x, where g{u) = J ^/a, and 
turn equation (|l.ip into the following system of equations 



/9 ON \ ut + ifiu) - y/a{u)q)x = bC[u], 

^ ' I q giu). = 0. 

Let us introduce the notation w = {u,q)' (here ' denotes the transpose), and write 



>.,w).M„,„.(M:,')^(«"'_;v;f5' 

Let us multiply each equation in (12. 8p by arbitrary Vu,Vq £ P^{Ii), integrate over 
the interval 7^, and use integration by parts, to arrive at 

dtuvu - I hu{w)dxVu + /i„(wj+i)i;^^^;^ - hu{wi)v^ ^ = b £[u]vu, 

li Jli ' "'/; 



h Jh 

where /i„(wi) = hu{ui,qi), Ui = u{xi), qi = qix^), v~^^ = Vu{x:[) and u^^ = Vuixf). 
Following Cockburn and Shu [7], we introduce the numerical flux 



(2.9) h(w-,w, ,-^ ^ ,_ .,+ 



^<?K 'O / V -5(ui) 



where n") = r/.c=( _^;;^ ^J^ 



1 /[F(«^ , _ 
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Ci2 = ci2(Wj^ , 'Wi ) is Lipschitz continuous in all its variables, and C12 = whenever 
a = or Wj^,w^ ~ 0. Note that cn > since / is an E-flux and, thus, the matrix 
C is semipositive definite. 

The LDG method consists of finding w = (u, g)', where 

fc k 

u{x,t) =^'^Up^i{t)ipp.i{x) and q{x,t) ^^^Qp^i{t)ipp^i{x) 
are functions satisfying 

(2.10) 'f ^'^ ' ' ■''^ 

qVq - I hq{u)d.j;Vq + hq{Ui+l)v^j^^ - hq{Ui)v+^ = 0, 

for all Vu, Vq e P'^{Ii), i e Z, and initial conditions for u and q given by (|2.6p . 

3. L^-STABILITY FOR NONLINEAR EQUATIONS 

We will show that in the semidiscrete case (no time discretization) both the DDG 
and LDG methods are L^-stable, for linear and nonlinear equations. 

In this section and the subsequent one, we assume the existence of solutions ii and 
w = (u, q)' of the DDG and LDG methods (|2.5p and (|2.10p . respectively, satisfying 
u,u,q G ^^([0, T]; V^ n L^(M)), in which case the integrals containing the nonlocal 
operator £[•] are aU weh defined. Indeed, by Lemma [Ql V^ n ^^(R) C iJ^/2(M), 
and hence all integrals of the form 



/ Lpi C[ip2] for y'l, <^2 e v'' n L^ 



can be interpreted as the pairing between tpi g H^/'^iS.) and C[ip2\ & H ^/^(M). 
Here i?"^/2(K) is the dual space of H^I'^iR), and C^] e H^^/'^{R) whenever 
ip e ff-^/2(R) (cf. CoroUary A.3 and proof in [5 ). 

Remark 3.1. The existence and uniqueness of solutions in C^{[0,T];V'' D i^(R)) 
can be proved using the Picard-Cauchy-Lipschitz theorem. The argument outlined 
in [SJ Section 3], can be adapted to the current setting since all numerical fluxes 
are (locally) Lipschitz (cf. [7] for the LDG case). For the DDG method with fc > 2, 
additional differentiability on a is needed for this proof to work. 

3.1. DDG method. Let us sum over all i € Z in (|2.5I) . integrate over t € (0,T), 
and introduce the functional 



MdDg[u,v]= / UtV~ y2\f('^i)['"i]+ / /(")«a; 

Jo JM Ja jg2 '- J Ii -' 

\h{ui)[v.i]+ I a{u)uxVx -hi / C[u]v. 
L Jii J Jo Jm 



(^■^^ rT . r . rT 



^^ez 



JO 

Let us define 



JO ._„ Jli Jo 






where 7 G (0, 1) and a > 0. Note that Tt > since a > and, using the Taylor's 
formula, [A(ui)][ui] — a(^i)[ui]^ > where and S,i £ [u{x~),u{x^)], i e Z. 
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Theorem 3.2. (Stability) Let u be a solution of (|2.5p such that both u,A{u) and 
their first k derivatives are sufficiently integrable. Then 

l|w(->r)|||2(K) + 2rT[M] + 6cA / |u(-,i)|^A/2(R)C?i< ||uo|li2(R). 

•^ 

Remark 3.3. Since u G C^{[0,T];V''nL'^{R.)) and /(O) = 0, all terms in ([321) below 
are well defined - except for 



I \y2h{ui)[ui]+ a{u){u,j:Y 

^ ,pz Jh 



When k > 2, additional integrability of u,A{u), and their first k derivatives, is 
required in order to give meaning to the /i-term. 

Proof. By construction, MddgIujv] = for all v E V'' H L'^{M.). If we set v ~ u, 
we obtain 

iitu- / ^l/(ui)[ui]+ / f{u)ux 
(3.2) 



iez 



+ / V \Hii^)[u^] + I a{u){uxf] -b f f £[u]u = 0. 
Jo ,,^2 ^ Jh J ^0 Jr 

Next, as a direct consequence of (j2.3p and a change of variables, we see that 



(3.3) 






f{u)Ua 



li 



<0. 



Since /i satisfies the expression (J2.4D . 

'■^ ^^^[4(_«, 

Ax 



iez 



^[ti]^ = -y|u|HA/2( 



E / a{u){uxf 



n 



(3.4) / E^^("»)N>«/ E^^N 

Finally, using Lemma lA. II 

(3.5) 

We conclude by inserting p.3p . p.4p . and p. 51) into 

3.2. LDG method. By summing over all i e Z, we can rewrite (I2.10p as 

/ dtUVu - E ^u{^i)[vu,i] + / hu{w)dxVu ) ^b £[u]vu, 

JTR jgg V Jli J JR 

/ 9^9 ^ E \k^i)\^qA + / hq{u)dxVq 1 = 0. 

We add the two equations and integrate over t e (0, T) to find Mlug[w, v] = for 

Mldg[w,v] = / / utVu+ / / qvq 
Jo Js. Jo Jm 

- / E fh(w,)'[v,] + / h(w)'a,v) -b [ [ C[u]v^, 

Jo fpr, \ Jli / Jo JR 



(3.6) 



where h(wj) = {hu{yfi),hq{ui))' , v = {vu,Vq)' and v^ = (««,», Wq,i)'- Moreover, let 
(remember that, as noted earlier, the matrix C is semipositive definite) 



•^0 JGZ 



DG FOR FRACTIONAL CONVECTION-DIFFUSION EQUATIONS 7 

Theorem 3.4. (Stability) // w = (u,g)' is a C^{[0,T];{V'' D L'^)'^) solution of 
([XTU)) . then 

•^ 

Here, as opposed to Thcorcni l3.21 no further integrability of the first k derivatives 
of the numerical solution w ~ (u, q)' is needed. The reason is that the numerical 
flux h has been built without the use of derivatives of w = (u, q)' . Each term in 
expression l\'3.7\i below is well defined thanks to p.8|) . the fact that /(O) = (which 
implies that cii(O) = 0), Ci2(0) =0, and w,ge C^{[0,T];V'' f] L^{R)). 

Proof. By construction, Mldg{^, v) = for aU v = (w„, u^)', Vu,Vq G V'' Ci L^{M.). 
We set V = w and find that 

cT I- i-T 



(3.7) 



/ [ utu+ f [ f~ [ y, f h(w»)'[w.] + / h(w)'9,w 
Jo JR Jo Jr Jo -^n, \ Jii 



^0 Jr 



Here we also used the fact that 

(3.8) - / ^ fh(w.)'[w.] + / h(w)'a,w') = / ^[w,]C[w,], 

Jo ^^^ \ Jh / ^0 i^j^ 

see [7] for a proof. To conclude, insert p. 81) and p. 51) into (|3.7p . D 

4. High-order convergence for linear equations 

In this section we consider the linear problem 

\ut + cUx=Uxx + 'bC[u] {x,t)£QT, 
]u{x,0) = uo{x) x G M, 

with the aim of proving that the DDG and LDG methods converge to a regular 
solution of (|4.ip with high-order accuracy. 

Lemma 4.1. Let uq E iJ'^+^(M). with k > Q. There exists a unique function 



€ H^+^Qt) solving gU, wh 



ere 



H^t'iQr) := {0 e L'{Qt) : ||9r9>llL2(Q,) < ^ for all < r + 2m <k + l}. 
Moreover, \\u{-,t)\\Hk+^R) < ||uo||h'=+i(k). 

Proof. Since the equation is linear, we can pass to the Fourier space. In view of 
(|1.2p . the Fourier transform of (14.11) is iit + i^cii — —^"^ii — h\£\^u. It follows that 

^(?,i)=^io(Oe"«^+«'^+''l«l'^ 

By the properties of the Fourier transform, the above expression implies the ex- 
istence of a unique L^-stable weak solution of (|4.ip . The L^-stability for higher 
derivatives can be obtained by iteration as follows: take the derivative of (|4.ip . use 
the Fourier transform to get stability, and iterate up to the fcth derivative. Regu- 
larity in time follows from the regularity in space since equation (|4.ip implies that 
d^u = {-cdx + dl + hCfu. n 

In the following two theorems we obtain L^-type error estimates for the DDG 
and LDG methods in the case that equation (|4.ip has iJp+^-regular solutions. 
(Note that the time regularity does not play any role here). To do so, we combine 
estimates for the local terms derived in [71 HO] with estimates for the nonlocal term 
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derived by the authors in [5 . In [6^ it was observed that most relevant numerical 
/ fluxes reduce to 

in the linear case. In this section we only consider this / flux. 

4.1. DDG method. 

Theorem 4.2. (Convergence) Let u G H^'^^^{Qt), k > 0, be a solution of (|4.1|) 
and u e C\[0, T];V'' D L'^{R)) be a solution of (ES]). With e = u-u, 

i-T 



bcxj^ |e|^v.(E) = OWAx^^ 



Remark 4.3. The error C'(l)Ax^'^ is due to the diffusion term Uxx- The errors from 
the convection term cu^ and the fractional diffusion term bC[u\ are of the form 
0(l)Aa;2'=+i and 0{l)Ax^''+^-^ respectively. 

Proof. Let us set 

Ma[u,v]= / utv+ / UxVx+ / y2h{ui)[v.i], 

Jo JR Jo Jr Jo jg2 

Mf[u,v] = ' f Y] \f{u,)[v,] + [ cuvX 

Mc[u,v] = -b / C[u]v. 

Jo Js. 

With this notation in hand, we can write p.ip as 

Mddg[u, v] = Ma[u, v] + Mf[u, v] + Mc[u, v]. 

Let Pe be the L^-projection of e into V'^, i.e., Pe is the y'^ni^(M) function satisfying 



(Fe{x) — e(x)'j(pji{x) dx — for all i E Z and j ~ {0, . . . , k}. 
It 

Note that Pe £ H^/'^{R) since ¥>" D L'^{R) C H^^'^{R) by Lemma El For all 
V eV'T] L2(M), we have Mddg[u, v] = since u is a DDG solution of (gH]), while 
MnDG[u,v] — since w is a continuous (by Sobolev imbedding) solution of (jl.ip 
and hence a solution of (|4.ip . Thus MDDG[e,v] = 0, and by bilinearity {h is linear 
since a = 1), 

(4.2) MoDGiPe, Pe] = MooGiPe - e, Pe]. 

One can proceed as in [20j (in that paper, combine the last inequality of the proof 
of Lemma 3.3 with Lemma 3.2 and (3.5)) to obtain 

(4.3) Ma[Pe-e,Pe] = i / / (Pe,)^ + i /" V /i(Pe,)[Pe,] + 0(l)Ax2'=. 

-^ Jo Jr ^ Jo jg2 

Moreover, proceeding as in O Lemma 2.17], 

(4.4) M/[Pe - e,Pe] = J^ / ^[Pe,]^ + Oil)Ax^''+\ 



DG FOR FRACTIONAL CONVECTION-DIFFUSION EQUATIONS 

As shown by the authors in [5], 

Mc\Pe - e, Pe] - Mc [Pe, Pe] = 6 / / £[e]Pe 

Jo Jr 

(4.5) =M fc[Fe]Ve+^f f C[e]e ^ ^ f f C[e -Ve]ie -¥e) 

^ jQ JwL ^ Jo Jm ^ Jo Jr 

.7 •'0 "^ 

where M£[Pe,Pe] = ^ J^ |Pe|^v2(R) fLemma [XT]) and 

(4.6) ||e-Pe|p^v.(M)<0(l)Aa;2fc+'"\ 
By (j3.ip . Lemma [A. 1[ and the definition of /, 



MDDG[Ve,Pe] = j{¥e\ + ^ / V[Pe,]2 + / ( (Fe^f 
Jr ^ Jo jg2 -^0 -^i^ 

+ / ^/i(Pe,)[Pe.] + ^/ |Pe|^v.(K)- 
Inserting this equation along with (I4.3p , (14.41) , and (14. 5p into (j4.2p then shows that 

m rri rri ^7^ 

/ /(Fe^)t + ^/ EP^^]' + |/ [i^^'^f + U EMPe.)[Pe.] 
and, using the admissibility condition (12.41) . 



^CA /■ 1^ ,2 , ^CA 



T 

4 j^ .--lffA/2(R)T-— y^ |e|^v2 



0(l)Ax 



2fe 



To conclude, we need to pass form Pe to e in the above expression. This has already 
been done for the diffusion term in Section 3 in |20j and for the convection term in 
the proof of Lemma 2.4 in [7]. For the nonlocal term, we see that by (j4.6p 

and the conclusion follows. D 

4.2. LDG method. 

Theorem 4.4. (Convergence) Let u E H!^'^^^{Qt), k > 0, be a solution of (|4.ip 
and w = {u,qY G C^([0,T]; V'' (1 i^) be a solution of ((^3)) . Mi/i e„ = u - u ond 

/ e^(z, T)+ f f el + eT[e] + 6ca / |e„|^v.(K) = 0(l)Aa;2'=. 
Jr Jo Jr Jo 

Proof. Let us choose a test function v — (u„,Wg)', Vu,Vq e V'^ n L^(R), and define 

M,[W,V] = / / Utt'u + / l^q- Y] ( h(w^)'[v.] + / h(w)'a^V 

Jo Jr Jo Jr Jo ^^^ ^ •^■f» ^ 
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With this notation at hand, we can write p.6p as 

Mldg [w, v] = Ml [w, v] + Ale [w, v] , 

where Mc is defined in the previous proof. Proceeding as in the proof of Theorem 
K2[ we find that 



(4.7) Mldg [Pe, Pe] = Mldg [Pe - e, Pe] . 
In [7] (Lemma 2.4) it is proved that 

(4.8) M,(Pe-e,Pe) = ieT[Pe] + i / f Pe^ + 0{l)Ax^'' 

By dSH), dXU, and Lemma Ell 

MLDGi^e, Pe] = / / (Pe2), + / / Pe^ + eT[Pe] + ^ / 
Jo Js. Ja JM ^ Jo 

By inserting this inequahty along with (j4.8p and (|4.5p into (14. 7L we find that 

fi^^Dt + lf /Pe2 + ieT[Pe] 
Jk ^ jo Jr ^ 

•^ •'0 



|2 



eulH^/2(R)- 



The conclusion now follows as in the proof of Theorem 14.21 D 

5. Convergence for nonlinear equations 

In the nonlinear case we will show that the DDG method converges towards an 
appropriately defined entropy solution of (jl.ip whenever piecewise constant ele- 
ments are used. In what follows we need the functions 

Vk{s) = \s-k\, 
Vkis) = sgn(s-fc), 
qk{s)^Vk{s)ifis)-f{k)), 
r,is) = 7j',{s){A{s)-A{k)). 

Remember that A{u) ~ J a, and let C^'^{Qt) denote the Holder space of bounded 
functions (f> : Qt — > R for which there is a constant c^ > such that 



\(j}{x,t)-(J3{y,T)\<c^^\x-y\ + ^/\i~^\^ for all {x,t),{y,T) € Qt- 
We now introduce the entropy formulation for (|l.ip . 

Definition 5.1. A function u G L°°{Qt) is a BV entropy solution of the initial 
value problem p.ip provided that the following conditions hold: 

(D.i) ueL\QT)nBV{QT); 

(D.2) A(u)€C1'5(Qt); 

(D.3) for all non-negative test functions ip G C^(M x [0,T)) and all fc G R, 

r]k{u)ipt + qk{u)ipx + rk(u)ipxx + 'q'f.iu) C[u]lp dx dt 

+ / •nk{uQ{x))Lp{Q,x)dx>Q. 
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This definition is a straightforward combination of the one of Wu and Yin [SD] 
(cf. also [13]) for degenerate convection-diffusion equations {b ~ 0) and the one of 
Cifani et al. [5] for fractional conservation laws (a = 0). By the regularity of (p and 
u and Lemma lA.li each term in the entropy inequality (D.3) is well defined. 



Remark 5.1. The L^-contraction property (uniqueness) for BV entropy solutions 
follows along the lines of [22], since the BV^- regularity of u and the L°°-bound on 



A{u)x makes it possible to recover from (D.3) the more precise entropy inequality 



utilized in [53] for L^ n L°° entropy solutions 

We will now prove, under some additional assumptions, that the explicit DDG 
method with piecewise constant elements (i.e., A; = 0) converges to the BV entropy 
solution of (jl.ip . In addition to convergence for the numerical method, this also 
gives the first existence result for entropy solutions of (II. ip . 

5.1. The explicit DDG method with piecewise constant elements. When 

piecewise constant elements are used (/c = in (j2.7p ). equation (j2.5p takes the form 



ut + fiui+i) - I{ui) - h{ui+i) + h{ui) = b C[u]. 
h J h 

Since u{x, t) — "Ylii^z Ui{t)\i{x) (i.e., (/jQ.i — li, the indicator function of the interval 
/i), we can and will use the admissible flux h{ui) = ■^^[^(ui)] (which satisfies (|2.4p 
with fc = and /?o = 1) to rewrite the above equation as 

For At > we set t„ = nAi for n = {0, . . . , iV}, T = tfq, and ^" = 4>{xi,tn) for any 
function (j). By a forward difference approximation in time, we obtain the explicit 
numerical method 

At Ax 

^^■^' AiUr+,)-AiUr) ^ Ajun - AjUr^,) _ b y^^.^,n 

Aa;2 Aa;2 Ax ^ ^ ^ ' 

where the weights G* = /^ -Cfl/J for all (i.,j) e Z x Z. All relevant properties of 
these weights are collected in Lemma [A. 21 Next we define 

D^V, = ±^ (C/,±i - U,) and £([/"), = ^ / C[U^] d^ = ^J2 '^)^h 



jei 



where [/" is the piecewise constant interpolant of [/": 

U'\x)^Ul\ xe[x^,x^+i). 
The explicit numerical method we study can then be written as 

f{ur, u^Xi) - D+A{ur)\ = bc{un., 



vr^-UT ^^_ 



(5.2) { ^* 

Uf^^JiM'^)dx. 

As we will see in what follows, the low-order difference method (j5.2p allows for a 
complete convergence analysis for general nonlinear equations of the form (jl.ip . 

Let us now prove that the difference scheme (|5.2p is conservative (P.l), monotone 
(P. 2), and translation invariant (P. 3). 
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(P.l) Assume C7" e L^{R) D BV{R). By LemmaEl] 



(5.3) ^^|qc/;| <^|/:[[7"(x)]|dx<c,q|L7"||^7('«)|c7"| 



l_By(R)i 



and hence we can revert the order of summation to obtain 

since J2iez^) = by Lemma (jA.21) . By summing over all i ^ Z on each 
side of (I5.2L we then find that 

(P. 2) We show that U"^^ is an increasing function of ah {[/"}igz. First note 
that 



dU. 



n+l 



— > for i^ j, 



since / is monotone and G* > for i 7^ j by Lemma IA.2I By Lemma IA.2I 
we also see that G* = —d\Ax^^''^ < 0, and hence 
a7j«+i At r - 

Here 9„. / denotes the derivative of /(ui, U2) w.r.t. u^ for i = 1, 2. Therefore 
the following CFL condition makes the explicit method (|5.2p monotone: 



(P. 3) Translation invariance {V^ ~ U'j^+i implies F/' = U^_^i) is straightforward 
since (|5.2p does not depend explicitly on a grid point a;^. 

Remark 5.2. For several well known numerical fluxes / (i.e. Godunov, Engquist- 
Osher, Lax-Friedrichs, etc.), we may replace 

\\9uif\\L°°{M) + \\du2f\\L°°(R) 

in the above CFL condition by the Lipschitz constant of the original flux /. 
In the following, we always assume that the CFL condition (15. 4p holds. 

5.2. Further properties of the expHcit DDG method (|5.2p . Define 

\MlHz)^Y.\^^\^ ||f/|Uoo(^)=sup|C/,|, and |t/|w(z) = JI I'^'+i " '^^l' 
i£Z *^^ tez 

Lemma 5.3. 

*j \\U''\\l^Z) < I|wo||li(r), 
tt) ||C/"|U~(z) < ||wo||l~(r), 
in) \U"\bv{Z) < \uo\BViR)- 

Proof. Since the numerical method (|5.2p is conservative monotone and translation 
invariant, the results due to Crandall- Tartar [51 [Tl] and Lucier [5S1 [T3] apply. D 



DG FOR FRACTIONAL CONVECTION-DIFFUSION EQUATIONS 



13 



For all (x, t) G i?f = [xi, x^-f i) x [i„, i„+i), let uax{x, t) be the time-space bilinear 
interpolation of C/", i.e. 



^A.(x,t)-c/r + (c^;;i-t^r) 



(5.5) 



X — llAX 



iAa 



iu, 



+ (t/^r - ur' - c/r+i + un 



X — iAx 
Ax 



t~nAt 

At 
t-nAt 

At 



Note that uax is continuous and a.e. difFerentiable on Qt- We need the above 
bilinear interpolation - rather than a piecewise constant one - to prove the Holder 



regularity in (D.2) We will show that the functions A{uax) enjoy Holder regularity 
as in (D.2)[ and then via an Ascoli-Arzela type of argument, so does the limit A{u). 
The following lemmas which are needed in the proof of Theorem l5.8[ are nonlocal 
generalizations of the ones proved in [T3] . In what follows we assume / S C^ (M) , 
and note that the general case follows by approximation as in [T3] . 

Lemma 5.4. 



(5.6) 






L-(Z) 



< 



f{U^,Ul,)-D+A{Uf)- 



A;— — oo jGZ 



L~(Z) 



(5.7) 



f{ur,ur+i)-D+A{un- E T.^P^ 



-oo jtz2 



BV(Z) 



< 






BV{Z) 



Proof. Inequality dSH). Let us start by defining V^ = ff EL-oo(C^fe " ^fc 
This sum is finite since [/" £ L^{Z) for all n > 0. If we use (j5.2p . we can write 

yn+l ^ 



(5.8) 



Here we have used that [/" e L^(Z) n BV{Z), f and A are Lipschitz continuous, 
and /(O) = to conclude that the sum EL-oo -^-[/(C/", t/"+i) - D+A{U^)] is 
finite and has value [f{Ul\ U^_^i) — D+A{Ul^)]. Next we rewrite the right-hand side 
of (US in terms of {F,"}^gz. By dEHD, 



-+i = vr [/(t/f, ^r+i) - f{ur\ui^+i') D+i^iun A{ur')) 



(5.9) 



k—~oo j^'L 



We prove that 



(5.10) Y. Y.G^m-ur') = ^T.G)vr- 

Indeed, note that D_V," = it (^" " C^""^) and 

E G'^i - E G^i^ = E ^r^' 

jez jez jGZ 
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since G*+i == G^ . Thus, 



E Y.G^(ur-ur')-^t y: eg^-k 



k— — oc ji 



k— — oo ji 



At 



£ E EG^(^-^i) 



Ax 

A^ E E^j'^"-A^ E E^j'^i-i 

/c— — oo jGZ /c— — oo jGZ 

^ EE^^^'-^ E E^r^ 



k— — C!C j'GZ 



i-1 



A2; '^-^ ^-^ -' -' At ^^— ^ ^-^ ^ ^ 
At 



k—~oo j?'GZ k— — oo j?'GZ 

Using Taylor expansions, we can replace the nonlinearities /, A with linear approx- 
imations as follows. We write 

(5.11) /(f/f , f/r+i) - f(ur\ c/r+i') = Ai/;",i?_c + At/£,i?_i-,^i, 

where />,, = 9i/«,C/f+i), Z^", = d2.f{Ur\al,,) and <,cif e {Ur^,U?)- 
Similarly, we write 



(5.12) 



A{u-) - A{ur') = a(A")(c^r - ur') = Ator^-^. 



where of = a(;3f) and /3f G (C/f'\C/f). Inserting (ICTIl and (I^1T|) - (IHT^ into 
expression (|5.9p returns 



(5.13) 



Ai 






or 



(5.14) 
where 



At 



yn+l ^ ^„^„ ^ ^ ^„^„ _^ ^„^„ ^ + ^ G; V^ 



iG2 



G" = 



^A + ^- 



At 

A^^ 



At 
Aa;2 



1-— (/M-/2".)-A^(«r + d) 

A^ .„ _At~ " 
Ax^'^'+i Ax ^2,. 



Since / is monotone and a > 0, A^, Gf > 0. Moreover, Bf + AlG^ > since the 
CFL condition (|5.4p holds true. Thus, since (|5.14p is conservative, monotone, and 
translation invariant (cf. the proof of Lemma [53| . ||V^"||l°°(z) < .•■ < ||^"^||l°°(z), 
and the conclusion follows from 
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Inequality dEZ]). Let us introduce Zf = V^" -V^-i- Note that, since G}l\ == G) 
for all (i,j) e Zx Z, 

iGZ jez jez 

Thus, (|5.13p can be rewritten as 



Ax 



jei 



or 



(5.15) 



Ai 



^„+i ^ ^„^„_^ ^ ^„^„ ^ ^„^„^^ + 7^ E ^5-^7' 



Ax 



iGZ 



where A", S", Cf have similar properties as A", B", Cf . Proceeding as in the first 
part of the proof, (|5.15l) can be shown to be conservative, monotone, and translation 
invariant. Thus \\Z^\\li(^z} < ■ ■ ■ < \\Z^\\l^(Z), and the conclusion follows from (j5.8p . 
We refer to [M] for the precise details concerning A^, -Bf , C". D 

The next lemma ensures that the numerical solutions are uniformly L-'^-Lipschitz 
in time (and hence BV in both space and time by Lemma [57 



Lemma 5.5. 

J2\ur~un< 



JSZ 



k— — oo j^l, 



At 

J\/(K) Ax 



\m — n\. 



Proof. Let us assume that m > n, the case m < n is analogous. Note that 



Eit^r-c^ri<EEi^^'-t^'i 



iGl 



l—n i^I 



m—l 



< 



A^EE 



l—n iGl 



£>_ 



f{Ui,UU,)-D+A{Ul) 



Ai 



■H^l 



jei 



Since D^ (EL-oo E,a G,^t^,) ^ i-^E.ezG^U, 



m- 



m—l 

^A^EE 

l—n zGZ 

To conclude, use (15.71). 



L>_ 



A:— — oo jGIa 



D 



We now show that the numerical solutions satisfy a discrete version of (D.2) 



Lemma 5.6. // |/(f/°, C/f+i) - D+A{U!) - EL-oo E,ez G,^C^,"|bv(z) < (X), i/ien 



\A{ur) - A(C/;)| = 0(1) [|i - j|Ax + Vk-^|Ai . 
Proof. Let us write 

|A(t/r) - ^(C/;)| < |^(C/r) - A{U^)\ + {AiUJ") - A{U^)\ =h+ h. 
We first estimate the term /i , then the term I2 ■ 
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Estimate of Ii. Using (I5.6p . (j5.3p . Leniina l5.3l ii). and the fact that / is Lipschitz 
continuous, 



\D+A{Un\ 



i ;ilL~(Z) 



< 






+ 



f{ur,u^i] 



L=°(Z) 



k— — oo jGZ 



L~(Z) 



L=°(Z) 



0(1). 



Hence h = 0{l)\i ~ j\Ax. 

Estimate 0//2. Take a test function (f) e C^(R), and let (f)i = (f){i/S.x). Let us 
assume m > n (the case tti < n is analogous). Using (|5.13p we find that 



Axj2q^uvr~v: 

l—n i^l. 

^^ E E <^^| (/m^-^/ + fLD-VU,) + D+{a-D^Vl) 



i& 



= Ax 



l—n i^ 



+ At 



EE^^E^jv^ 

i=n iez jez 



C1+C2. 



We use summation by parts to move D_|_ onto (j)i and the fact that fl\ /^^^ a" and 
\V^\bv(T) are uniformly bounded to arrive that 

Ci = 0{l)At{m - n) (||0|U=.(m) + W\\l^(^)) ■ 
For more details, see [H]. Then by (jO]), E»ezEjez l^'j^j I = ^(1)' and hence 

771—1 

(5.16) C2 < At||0iUo.(R) E EE i^^-^'i ^ 0(i)At("^ - ri)ii0iUo=(K). 

Therefore, 



-n 2^Z jG2 



(5.17) 



AxY,^dVr-V[' 



i& 



0{l)At{m-n) I 



The above inequality is exactly expression (40) in |14j . From now on the proof 
continues as in [14]. Loosely speaking we take an appropriate sequence of test 
functions (f>^ e C^ (M) to deduce from (|5.17p that 

AxJ2\Vr - Vn = OilW{m^n)At. 



By (j5.8p . Lemma [5.51 and inequality (j5.16p we also find that 

Ax J2 \Vr - CI = 0{l){m - n)At + AxJ2\D+A{ur} ~ D^A{U:\ 



iei- 



iei 
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and hence AoiX^iez l-D+^(f/") - D+A{U:^)\ = 0(l)^(m - n)Ai. We conclude by 
noting that 



/2 = |A(t/f)-A(C/")HAa; 



£ D+A{Ur)- J2 ^+^(^" 



< Ax^ \D+Aiun - D+Aiun\ = OilWim-n)At. 



iez 



Next we show that the numerical method 



n 



satisfies a cell entropy inequality, 



which is a discrete version of (D.3) 

Lemma 5.7. LetkeM. and ryf = |C/f - k\. Then 

(5.18) v7^' - V? + AtD^Q^ - AtD^D+\A{un - A(fc)| < Atr,',{Ul^+')C{U^^),, 

where Qf = /(C//' V fc, C/f+i V k) - f{Ul' A fc, L/f+i A k). 

Proof. Let us introduce the notation a/\h = min{a, fe} and aM b — max{a, &}. Note 
that ?7" = (t/" V fc) — (C/" A A;). Since the numerical method (|5.2I) is monotone, 

(£/r+^ V fc) - ({/f V k) f{Ur V fc, ;7r+i V fc) - /({/;Li V fc, L/f V fc) 
At Ax 

A{U^^^ V fc) - A(L/;' V fc) ^(C/f V fc) - v4(L/;Li V fc) 



Aa;2 



Aa;2 
<Ail(,,+^)(f/f+i)£(C/"), 



and 



{U^+^ A fc) - (£/f A fc) ^ /(J7r A fc, t/f+i A fc) - /(J7f_i A fc, t/f A fc) 
At Ax 

A{U^_^^ A fc) - A(L/f A fc) A{U'^ A fc) - A(C/f_i A fc) 



Aa;2 



To conclude, subtract the above inequalities. 



Aa:2 



n 



5.3. Convergence of the DDG method. We are now in position to prove con- 
vergence of the fully explicit numerical method (j5.2p to a BV entropy solution of 
(II. ip . Let us introduce B (cf. jH]), the space of all functions 2; : M ^- R such that 



f{z) - d^A{z) 



C[z] 



< 00. 



BV( 



In the following theorem we choose the initial datum to be in L^(R) H BV{M.) D B, 
which is done to make sense to the right-hand side of (|5.7p . Note that whenever 
z e L^R) n BV{R.), C[z] e L'^iR) by LemmaEH and hence 

d r 



C[z] 



BV( 



dx 



\C[z] 



< 00. 



Theorem 5.8 (Convergence for DDG). Suppose uq e L^{M.) n BV{R) n B, and 
let iiAx be the interpolant (j5.5p of the solution of the explicit DGG scheme (|5.2I) . 
Then there is a subsequence of {uax} O'nd a function u € L^{Qt) H BV{Qt) such 
that (a) U/\x -^ u in L\^^[Qt) as Ax — ^ 0; (b) u is a BV entropy solution of (|1.1D . 



Corollary 5.9 (Existence). If uq e ^^(R) n 

entropy solution of (jl.ip . 



n B, then there exists a BV 
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Proof of Theorem \5.8[ We will prove strong Lj^^ compactness, and hence we need 
the following estimates uniformly in Ax > 0: 

i) \\uax\\l-^{Qt) < C", 
n) WuAxWsviQT) < d- 
Estimate i) is a consequence of Lemma l5.3l and (j5.5p . while estimate ii) comes from 
the following computations (cf. _14_ for more details) . Using the interpolation ([57 
we find that 



N-l 






N-l 



Y.Y.Pr-,\'-ur'\ 



n=0 i& 



< T\U \bv(i)- 
Note that Lemma 15.31 has been used in the second inequality. Similarly, 

/ i"*i ^ ^ E E pr' u-\ + f^ E E l^m - u^ 



<r 



/(c7°,t/,%)-D+A(t/°)- Y. T.Gpj 



BV(Z) 



h— — oc j(zZ 

where Lemma 15.51 has been used in the second inequality. Hence, there exists a 
sequence {uAxi}i£N which converges in L\^^{Qt) to a limit 

ueL^{QT)f^BV{QT). 



Next we check that the limit u satisfies (D.2) We define wax = A{uax)- Note 
that A{uAx) -^ -^iu) a.e. since uax -^ u a.e. (up to a subsequence) and A is 
continuous. Now choose (a;, t), {y, r), (j, n), (i, m) such that (a;, t) G -R" and (y, r) € 
i?™ for Rf = [xi,Xi+i) X [t„,t„+i). Then, 

\wAxiy, t) - WAx{x, t)\ < \wAx{y, t) - WAx{iAx, ■mAt)\ 

+ \wAx{iAx,mAt) — WAx{jAx,nAt)\ 

+ \wAxijAx,nAt) -WAx{x,t)\ 

= h +h+h- 



Note that by Lemma [5T6l I2 — 0(l)(|i — j|Aa;+ y'\m — n|At), while by Lemma [576 
again, ([53]), and A' ^ a e L°° , h + I3 = 0{l){Ax + y/At). Thus 



\wAxiy,T)-WAxix,t)\ =0(1) \\y -x\ + Vk -t\ + Ax + VAi 

We also have that wax = A{uax) is uniformly bounded since A is Lipschitz and 
iiAx is uniformly bounded. By essentially repeating the proof of the Ascoli-Arzela 
compactness theorem, we can now deduce the existence of a subsequence {wax} 
converging locally uniformly towards the limit A{u). By the estimates on wax, it 
then follows that 



(5.19) 



Aiu)eC'^HQT). 



Finally, let us check that the limit u satisfies (D.3) in Definition 15.11 Here we 
need to introduce a piecewise constant inteporlation of our data points t/". We call 

UAx{x,t) ^ Ul for aU (a;,i) G [xi,Xi+i) x [t„,i„+i). 

We do this since the discontinuous sign function ry^ makes it difficult to work with 
the bilinear interpolant uax in what follows. The need for the piecewise linear in- 
terpolation was dictated by the condition (D.2) continuity of the functions A(uax) 
were needed to prove Holder space-time regularity for the limit A[u) (cf. the proof of 
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(|5.19p ). To verify that the limit u also satisfies (D.3) the piecewise constant inter- 
polation uak suffices since, as we already have strong convergence for the piecewise 
linear interpolation, strong convergence toward the same limit u for the piecewise 
constant interpolation is ensured thanks to the fact that 

We now take a positive test function (^ G C^(IR. x [0,T)), and let (^" = (p{xi,tn). 
We multiply both sides of (jS.lSp by ff, and sum over all {i,n). Using summation 
by parts, we obtain 

N-l „+l 

n=0 ieZ 



A-AtEE^r^^-^p^ 



JV-l 
JV-1 

(5-20) + AxAt J2 E l^(^") ~ Mk)\D-{D+^:) 

JV-1 

+ Ax ^ ^077.0 >0. 

zSZ 

A standard argument shows that all the local terms in the above expression converge 
to the ones appearing in the entropy inequality (D.3)[ see e.g. [THKIJ. Let us look 
the term containing the nonlocal operator £(■). We can rewrite it as 
rT+At j- 

I r]'k{'>J-Ax)C[uAx]<f'dxdt + R, 

At JR 

where R — > and ip is the piecewise constant interpolant of 99". Indeed, let us 
write 77^([/f+i)£(C/"), = 7?^,(L/f+i)/:(L/" - L/"+i), + ry^(L/f+i)/:(C/"+i>,. Note that 

Af-l Af-1 „ 

AxAi^^ !£([/" -[/"+i),|^r<Ai||^|Uo.(Q,)^ / \C[U^\x)-U"+\x)]\dx, 

n=0 iGZ n=0 "^^ 

where the last quantity vanishes as Ax — ^ by L^-Lipschitz continuity in time (cf. 
Lemma [5T51 and also Lemmas 15.31 and IA.1[) . Next, 

N 

N N 

n=0 iGZ n=0 iGZ 

where the first term on the right-hand side vanishes as Ax — >■ since there exists a 
constant c^ > such that |(/3" — V'""''^! — Ci^Aa; for all (i, n). To conclude, we prove 
that up to a subsequence and for a.e. /c S M, 

rT+At (. /. 

(5.21) / / rj'i.{uAx)C[uAx]'^dxdt — > / irj'j^{u)C[u]ipdxdt. 

J At Jr Jqt 

This is a consequence of the dominated convergence theorem since the left hand 
side integrand converges pointwise a.e. to the right hand side integrand. Indeed, 
first note that ip ^ ip pointwise on Qt, while a.e. up to a subsequence, u^x —>■ u 
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on Qt- We also have ?7^('Sax) ^ "Hki"^) ^-^^ i'^ Qt since for a.e. fc S M the measure 
of {{x,t) G Qt ■ u{x,t) — k} is zero and 77^ is continuous on M\{A;}. FinaUy if 
the (compact) support of ip is containd in [—B.,R] x [0,r], i? > 0, then a trivial 
extension of Lemma lA.ll implies that 



/ \£[uAx - u]\ dx dt < C),C \\uAx - u\\\i^_„ jjAuAx 

J\-R.B]x\0,T] Jo y ' t 



ll-A L-,. ,,|A 

'R.R]xlO,T] '- " " ■' 



where the last quantity vanishes as Aa; — >■ since uax — >■ u in L\oc{Qt)- Then 
C-[u\x] ^ i^[u] a.e. in [—R,R] x (0,T) up to a subsequence. Convergence for all 
fc e R can be proved along the lines of [2S1 Lemmas 4.3 and 4.4]. D 

5.4. Remarks on the LDG method. The derivation of the LDG method in the 
piecewise constant case is not as straightforward as the one for the DDG method. 
Indeed, the numerical fluxes introduced in (j2.9p depend on the choice of the function 
C12, and computations cannot be performed until this function has been defined. 
Our aim now is to show that the LDG method reduces to a numerical method 
similar to (|5.1I) for a suitable choice of the function C12 . 

Let us for the time being ignore the nonlinear convection and fractional diffusion 
terms and focus on the problem 



ut ~ dx\/a{u)q = 0, 
q - dxg{u) = 0, 
u{x, 0) = uq{x). 

The LDG method ((2T0l) then takes the form 

{/, Ut + /i„(wi+i) - ^„(wj) = 0, 
JU ^ + hq{Ui+l) - hq{Ui) = 0, 

where u{x, t) = J2i(iZ Ui{t)lii {x), q{x, t) = Y.i& Qi{i)'^h (a^ ),and the fluxes (/i„, hq) 
are defined in (j2.9p . Let us insert ii and q into the system (j5.22p . and use the flux 
(EH) to get 

r |C/.Ax - Mi^titim Q.±^ _ ci2(Q.+i - Q.) 



(5.23) 






' 2 

Let us choose the function C12 to be 

(5.24) ci2(t/.,C/.-i)^- ^^_^^_^ . 

Inserting (|OI)) into (|05)) then leads to 

dt^«^-^ Ui+i-Ui Vj+1 ^ u.-u._-^ Vj — o, 

^ g(J7,)-g(J7.-i) 
Ax ' 



or 



dt ' Ax C/,+1 - [/, Aa; t/, - C/,_i 
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For the full equation (|l.ip , this choice of C12 along with a forward difference ap- 
proximation in time, lead to the following piecewise constant LDG approximation: 

ur+'-ur f{u?.u?+i)'f{u?-i^ur) 

At Ax 

Aa;2 up.. - up Ax^ UP -UP , Ax ^ ^ ^ ' 

Remark 5.10. We do not prove convergence for the numerical method (j5.25p . How- 
ever, we note that 

dgV dA 



^du J du^ 
since g = J \/cl and A = J a. Roughly speaking this means that 

{gm+i) - 9{up)r 



U,+i - u. 



A{UP^,)-A{UP), 



and hence that (|5.ip and (j5.25p are closely related. Experiments indicates that the 
two methods produce similar solutions (cf. Figure [3]). 



6. Numerical experiments 

We conclude this paper by presenting some experimental results obtained using 
the fully explicit (piecewise constant) numerical methods (|5.2p and (|5.25p . and the 
DDG method (|2.5p with fully explicit third order Runge-Kutta time discretization 
and piecewise constant, linear, and quadratic elements. In the computations we 
have imposed a zero Dirichlet boundary condition on the whole exterior domain 
{|x| > 1}. In all the plots, the dotted line represents the initial datum while the 
solid one (or the dashed-dotted one in Figure |31) the numerical solution at i = T. 

Remark 6.1. The operator C[u\ requires the evaluation of the discrete solution u on 
the whole real axis, thus making necessary the use of some localization procedure. 
In our numerical experiments we have confined the nonlocal operator C[-] to the 
domain Vl — {\x\ < 1}. That is to say, for each grid point {xi,tn) G il x (0,T) we 
have computed the value of it at time tn+i by using only the values u{xi,tn) with 
Xi G fl. 

We consider two different sets of data taken from [M" . In Example 1 we take 



ai{u) 



(Ex.1) 



"0,1(2;) 



«^ 




{' 


for u < 0.5 


<2.5u-1.25 for 0.5 <u< 0.6 


[0.25 


for u > 0.6, 


'0 


for X < -0.5 


5x + 2.5 


for - 0.5 < X < -0.3 


1 


for - 0.3 < X < 0.3 


2.5 -5a; 


for 0.3 < a; < 0.5 





for X > 0.5. 
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In Example 2 we choose 

/2 = ^/l, 




0-2 = 4ai, 
(Ex.2) /^ 

"0,2(2;) = < -2.5a; 



for X < -0.4 

for - 0.4 < a; < 

for X > 0. 


Furthermore, in Example 3 we use 
.^„o^ a,{u) = 0.1, 





"0,3(2;) = exp f- (^—j y 



The numerical results are presented in Figure [TJ [21 El and 21 The results confirm 
what we expected: the solutions of the initial value problem (jl.ip can develop 
shocks in finite time (this feature has been proved in [3] for the case a = 0). In 
Figure [1] and JH you can see how the presence of the fractional diffusion C influences 
the shock's size and speed. In Figure JH you can see how the accuracy of DDG 
method (|2.5p improves when high-order polynomials are used (k = 0, 1, 2). 





(a) Ut + f{u)^ = {a{u)Uj:)a: 



(b) Equation 11. Il l with A = 0.5 



Figure 1. ((E^ : T = 0.15 and Ax = 1/640. 

In Figure |3l the dashed-dotted curve represents method (j5.2p . while the solid one 
represents method (|5.25p . The two numerical solutions stay close, and numerical 
convergence has been observed for finer grids. Note that here we have set 6 = (no 
fractional diffusion) in order to stress the differences between the two methods. 

The numerical rate of convergence for the solutions in Figure [H (b), [51 (b), andJH 
(b) are presented in Table [H We have measured the L^-error 

Eax.p = \\uAx{-,T) — Ue{-,T)\\j^p^^y 

where Ug is the numerical solution which has been computed using a very fine grid 
(Aa; = 1/640), the relative error 



Rax,p 



1 



ll"e(-,T)||^,( 



Eax,p, 



and the approximate rate of convergence 

aAx..p^ i-^ ^j (l0g^Ax,p-l0g^Ax/2,p 
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(a) ut + f{u)x = {ci{u)ux)x (b) Equation 111. Il l with A = 0.5 

Figure 2. (|E^ : T = 0.25 and Aa; = 1/640. 



(a) T = 0.0625 



(b) T = 1 



Figure 3. (jEx.2p : solutions oiut + f{u)x ~ ia{u)ux)x at different 
times using methods dOJ and (|5T25)) (Ax = 1/160). 



Table 1. Error, relative error, and numerical rate of convergence 
for the solutions in Figure[T] (b),[2](b), and|l](b). 





Fi 


gure [11(b) 


Fi 


gure[l(b 


) 


Fi 


gure 11(b) 




Ax 


Eax,1 


^Ax.l 


"Ax,! 


Eax.I 


RAx,1 


OiAx.l 


Eax.2 


RAx.2 


aAx.2 


1/10 


0.0706 


0.0942 


0.97 


0.0474 


0.0550 


0.86 


0.009000 


0.093595 


2.00 


1/20 


0.0361 


0.0482 


0.92 


0.0261 


0.0302 


0.49 


0.002300 


0.023493 


1.85 


1/40 


0.0191 


0.0255 


0.57 


0.0186 


0.0216 


0.52 


0.000626 


0.006518 


1.54 


1/80 


0.0128 


0.0171 


0.60 


0.0130 


0.0150 


0.42 


0.000216 


0.002248 


1.10 


1/160 


0.0084 


0.0113 


0.76 


0.0097 


0.0112 


0.77 


0.000101 


0.001052 


1.04 


1/320 


0.0050 


0.0066 


- 


0.0057 


0.0066 


- 


0.000049 


0.000510 


- 



Our simulations seem to indicate numerical convergence of order less than one for 
the solutions depicted in Figure [1] (b) and [H (b) (nonlinear equations and piecewise 
constant elements), and numerical convergence of order higher than one for the 
solution depicted in Figure H (b) (linear equation and piecewise linear elements). 
In the last case we do not seem to reach the expected value 2 (cf. the statement of 
Theorem 14.21) . This deterioration of the numerical order of convergence for high- 
order polynomials has already been observed by the authors in [5]. The reasons 
behind this deterioration are still not clear. 
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(a) Piecewise constant {k = 0) with Ax = 1/20 (b) Piecewise linear {k = 1) with Ax = 1/20 





l\\ 

1 \ \^ 




(c) Piecewise quadratic (k = 2) with Ax = (d) Solution computed using Ax = 1/640 
1/20 

Figure 4. ([Ex^j) : solutions at T = 0.1 using fc = 0, 1, 2. 



Finally, let us remind the reader that no general results concerning the rate 
of convergence of numerical methods for nonlinear equations like (jl.l[) have been 
produced so far. For more details, cf. [3]. 

Appendix A. Technical lemmas 

In this appendix we state some technical results from 5, that are needed in this 
paper. All proofs can be found in ISl. 

Lemma A.l. Let ip,(j) ^ L^(M) n BViJi). Then there exists C > such that 



(A.2) 



(A.3) 



<l>C[ip] 



ipC[ip\ = - 



^£[0], 



c\ 



\z — x\ 



1+A 



dz dx. 



Moreover, the last two identities also hold for all functions (f>,(p E H^''^{M.). 



To prove inequality (jA.l|) one can split the nonlocal operator £[•], using an auxil- 
iary parameter e > 0, into the sum of £e[-], the operator containing the singularity, 
and C'^[-], the remaining part of the original operator. The operator Ce[-] can then 
be treated using the control on the bounded variation, while the control on the L^- 
norm is needed for the operator £^[-]. To obtain exactly estimate (jA.ip the optimal 
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value of e must be chosen. The proof of (|A.2p - and thus of (jA.3|) - is essentiahy a 
change of variables. 

Lemma A. 2. For all ii,i) G Z x Z, 

^ |G^| < oo, ^ G^ = 0, G) = Gi and G%\ = G). 

keZ feGZ 

Moreover, G' > whenever i ^ j , while 

(A.4) Gl^-c,\ I -^+ I -^ ) Ax'-' < 0. 




2 >1 



Lemma Fa. 2 1 is essentially a consequence of the form of the operator £[■] itself, 
and properties (jA.ip and (jA.21) . Property (IA.4I) comes from a precise evaluation of 
the integral G-. 



Lemma A.3. If<l)eV''n i^(M), t/ien G i/-(M) /or aZ^ A G (0, 1), and 

^Ht(R) - Ax' 



(A-5) ll'^lllA_<^ll'/'lli.( 



Lemma [A.SI is essentially a consequence of the fact that is a piecewise polyno- 
mial. The control on the i^-norm together with the piecewise structure of 4> ensure 
that its quadratic variation is bounded. Then, the finite quadratic variation plus 
the fact that (p is differentiable inside each interval li return (|A.5p . 
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